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ABSTRACT 

We present a sample of 8498 quasars with both SDSS ugriz optical and UKIDSS 
YJHK near-IR photometric data. With this sample, we obtain the median colour-z 
relations based on 7400 quasars with magnitude uncertainties less than O.lmag in all 
bands. By analyzing the quasar colours, we propose an empirical criterion in the Y — K 
vs. g — z colour-colour diagram to separate stars and quasars with redshift z < 4, and 
two other criteria for selecting high redshift quasars. Using the SDSS-UKIDSS colour- 
z relations, we estimate the photometric redshifts of 8498 SDSS-UKIDSS quasars, 
and find that 85.0% of them are consistent with the spectroscopic redshifts within 
|Az| < 0.2, which leads to a significant increase of the photometric redshift accuracy 
than that based on the SDSS colour-z relations only. As two tests, we compare our 
colour selection criterion with a small UKIDSS/EDR quasar/star sample and a sample 
of 4671 variable sources in the SDSS Stripe 82 region with both SDSS and UKIDSS 
data. We find that they can be clearly divided into two classes (quasars and stars) 
by our criterion in the Y ~ K vs. g — z plot. We select 3834 quasar candidates from 
the variable sources with g < 20.5 in Stripe 82, 826 of them being SDSS quasars and 
the rest without SDSS spectroscopy. We estimate the photometric redshifts for 3519 
quasar candidates with all UKIDSS YJHK data and find an accuracy of 87.5% within 
|Az| < 0.2 with the spectroscopic redshifts of 819 SDSS-UKIDSS identified quasars 
among them. We demonstrate that even at the same spectroscopy limit as SDSS, with 
our criterion we can at least partially recover the missing quasars with z ^ 2.7 in 
SDSS. The SDSS identified quasars only take a small fraction (21.5%) of our quasar 
candidates selected from the variable sources in Stripe 82, indicating that a deeper 
spectroscopy is very promising in producing a much larger sample of quasars than 
SDSS. The implications of our current results to the future Chinese LAMOST quasar 
survey are also discussed. 
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1 INTRODUCTION 

The number of quasars has increased substantially in the 
last decade primarily due to two large optical sky surveys, 
namely the Two Degree Field (2DF) survey (Boyle et al. 
2000) and the Sloan Digital Sky Survey (SDSS) (York et al. 
2000). 2DF has obtained spectroscopy of more than 20,000 
quasars (Groom et al. 2004), and SDSS has identified more 
than 100,000 quasars (Schneider et al. 2007; Abazajian et 
al. 2009). Quasar candidates in these surveys were mainly 
selected based on their optical colours derived from the pho- 
tometric data. 2DF mainly selected lower redshift {z < 2.2) 
quasars with UV-excess {u — hj <-0. 36) (Smith et al. 2005), 
while SDSS adopted a multi-band optical colour selection 
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method for quasar selection mainly by excluding the point 
sources within the stellar locus of the colour-colour diagrams 
(Richards et al. 2002). 90% of SDSS quasars have lower red- 
shift [z < 2.3), although some dedicated methods were also 
developed for discovering high redshift quasars (Fan et al. 
2001a, b; Richards et al. 2002). The lower efficiency in iden- 
tifying quasars with redshift between 2 and 3 is obvious in 
SDSS (Schneider et al. 2007), because these quasars usually 
have similar optical colours as stars and thus are mostly ex- 
cluded by the SDSS quasar candidate selection algorithm 
(Warren et al. 2000). 

In order to obtain a more complete sample of quasars, 
obviously we need to improve the previous quasar candidate 
selection methods. An important way to do this has been 
suggested by using the infrared K-band excess to identify 
the 'missing' quasars with redshift around 2.7 based on the 
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UKIRT (UK Infrared Telescope) Infrared Deep Sky Survey 
(UKIDSS) (Warren et al. 2000; Hewett et al. 2006; Maddox 
et al. 2008) □ Although the z ~2.7 quasars have similar op- 
tical colours as stars, they are more luminous in the infrared 
K-band. In addition, combining the optical colours in SDSS 
with the infrared colours in UKIDSS, it should be more ef- 
ficient to separate stars from both lower redshift quasars 
(z <3)(Chiu et al. 2007) and higher redshift ones (z >6) 
(Hewett et al. 2006) in the colour-colour diagrams. This ad- 
vantage needs to be confirmed by future spectroscopic ob- 
servations on the larger sample of quasar candidates selected 
with the SDSS and UKIDSS photometric colours. 

Multi-band colours are also crucial for the photometric 
redshift estimations of quasars. Based on the SDSS colour-z 
relations, about 70% of the photometric redshifts of quasars 
are consistent with the spectroscopic ones with the differ- 
ence |A2| < 0.2 (Richards et al. 2001; Wu et al. 2004; We- 
instein et al. 2004). This percentage can increase to 85% 
or higher if the extra GALEX UV photometric data (Ball 
et al. 2007) or Spitzer IRAC data (Richards et al. 2009b) 
are added. It has been demonstrated that the photo-z accu- 
racy should be improved if combining the UKIDSS infrared 
photometry with the SDSS optical photometry for galaxies 
(Maddox et al. 2008). Such improvement can be also ex- 
pected for quasar candidates with both SDSS and UKIDSS 
photometric data, but has not been confirmed by the photo- 
metric redshift estimations on a large SDSS-UKIDSS quasar 
sample so far. In addition, the reliable photometric redshift 
estimations based on multi-band colours are also important 
for preparing the quasar candidates for the future spectro- 
scopic surveys (Richards et al. 2009a). 

Large quasar survey is one of the key projects of the 
Chinese Large Sky Area Multi-Object Fibre Spectroscopic 
Telescope (LAMOST), which is a novel reflecting Schmidt 
telescope with 4 meter effective mirror size, 20 square degree 
field of view (FOV) and 4000 fibres (Su et al. 1998). As the 
most efficient optical spectroscopic telescope in the world, 
LAMOST finished its main construction in 2008 and has 
entered the commissioning phase. A pilot survey and the 
regular survey have been planned in 2010 and 2011-2015 
respectively. Unlike SDSS, LAMOST does not have its own 
photometric survey. Therefore, the input catalogue of LAM- 
OST quasar survey will largely rely on the photometric data 
from other existing surveys. Because of the large overlap of 
the surveyed area between SDSS and UKIDSS, a combina- 
tion of SDSS and UKIDSS photometric data is expected to 
help us to efficiently select a large catalogue of quasar can- 
didates and provide reliable photometric redshifts for the 
LAMOST quasar survey. 

In this study, we first present a large quasar sample with 
both SDSS and UKIDSS data, and then investigate whether 
there is a more efficient selection criterion for quasars and 
whether the combination of UKIDSS data can improve the 
photometric redshift estimations. As two tests, the results 
are applied to a small UKIDSS/EDR quasar/star sample 



1 The UKIDSS project is defined in Lawrence et al. (2007). 
UKIDSS uses the UKIRT Wide Field Camera (WFCAM; Casali 
et al. 2007) and a photometric system described in Hewett et al. 
2006. The pipeline processing and science archive are described 
in Hambly et al. (2008). 
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Figure 1. Upper panel: Histogram of the position offsets be- 
tween SDSS and UKIDSS for 8498 quasars. Middle panel: His- 
togram of the quasar redshifts. Lower panel: Histograms of the 
redshifts of quasars with different UKIDSS mergedclass parame- 
ter. Solid, dotted, long dashed and short dashed histograms cor- 
respond to mergedclass=-l (star), -2 (probable star), +1 (galaxy) 
and -3 (probable galaxy), respectively. 

and a sample of the variable sources in the SDSS stripe 82 
region. The implications of our study to the future LAMOST 
quasar survey are also discussed. 



2 A SAMPLE OF QUASARS WITH BOTH 
SDSS AND UKIDSS DATA 

We cross-identify all quasars in SDSS Data Release 7 (DR7) 
with the UKIDSS Data Release 3 (DR3), by finding the 
closest counterparts within 3" between the positions in two 
surveys and requesting all detections in both SDSS ugriz 
and UKIDSS YJHK bands for each quasar. To do the cross- 
identifications, we use the CrossID form available at the 
UKIDSS WFCAM Science Archive websitfl and use only 
the data in UKIDSS Large Area Survey (LAS) in order to 
avoid the mis-identifications in the crowded fields with lower 
Galactic latitudes. This results in a sample of 8498 quasars 
with both SDSS and UKIDSS data. 

In the upper panel of Fig.l, we show the distribution 
of the position offsets between SDSS and UKIDSS for 8498 



^ http:/ /surveys. roe. ac.uk:8080/wsa/crossID_form.jsp/ 
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Table 1. Parameters of 8498 SDSS-UKIDSS quasars 
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Note: All magnitudes arc in Vega system. The full table of 8498 quasars is available in the electronic version of the paper. 



quasars. Although we use the search radius 3" for identi- 
fying the UKIDSS counterparts of SDSS quasars, we find 
that 99.6% of the closest counterparts are within 0.5" of the 
SDSS positions. Only 33 identifications have offsets larger 
than 0.5", and only three of them have the offset larger than 
1" . Therefore, considering the subarcsecond position accu- 
racy in both SDSS and UKfDSS photometry, the chance 
probability of misidentification should be very small. The 
middle panel of Fig.l shows the redshift distribution of 8498 
quasars. The redshifts are distributed in a range from O.Of 
to 5.29, with 7468 quasars (87.9% of the total sample) hav- 
ing 2 <2.3 and 1030 quasars (12.1%) having z >2.3. Only 
five quasars have redshifts larger than 5. This redshift dis- 
tribution is quite similar as that of the whole quasar sample 
in SDSS (Schneider et al. 2007). The deficiency of z ~ 2.7 
quasars can be also clearly observed. In the lower panel of 
Fig. 1 we show the redshift distribution of the quasars with 
different UKIDSS mergedclass parameter. Among the 8498 
SDSS-UKIDSS quasars, 4905, 1181, 2383 and 29 quasars are 
classified with UKIDSS mergedclass=-l (star), -2 (probable 
star), +1 (galaxy) and -3 (probable galaxy), respectively. 
It is clear that most quasars with z < 0.5 are classified as 
extended in UKIDSS morphology, which confirms the signif- 
icant contributions of host galaxy light to the quasar emis- 
sions in near-IR (Maddox & Hewett 2006; Chiu et al. 2007). 
However, the fraction of quasars classified as point sources 
in UKIDSS morphology becomes dominant at 2 > 1. This is 
understandable because it is hard to detect the host galaxies 
of quasars at high redshift even in the near-IR band. 



In Fig.2 and Fig.3 we show the distributions of magni- 
tudes and magnitude uncertainties of 8498 quasars in both 
SDSS ugriz bands and UKIDSS YJHK bands. AU magni- 
tudes are in Vega system throughout this paper if not spec- 
ified. We convert the SDSS AB magnitudes to Vega magni- 
tudes using the following scalings (Hewett et al. 2006): u = 
u{AB) - 0.927, g = giAB) + 0.103, r = r{AB) - 0.146, i = 
i{AB) - 0.366 and z = z{AB) - 0.533. The magnitudes in 
SDSS are corrected for the Galactic extinctions using the 
extinction map of Schlegel et al. (1998). The quasar mag- 
nitudes are mostly in a range from 17 to 20 in all SDSS 
optical bands and UKIDSS Y and J bands, from 16 to 
19 in H band and from 15 to 18 in K band. The median 
magnitudes are 18.3,19.11,18.70,18.37, 18.13 for ugriz and 
18.16,17.79,17.22, 16.59 for YJHK, respectively. The mag- 
nitude uncertainties are typically smaller than 0.05mag in 
SDSS gri band, O.lmag in u band, and 0.15mag in SDSS 
z and UKIDSS YJHK bands. The median values of mag- 
nitude uncertainties are 0.03,0.01,0.01,0.01, 0.04 for ugriz 
and 0.03,0.03,0.04, 0.04 for YJHK, respectively. The me- 
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Figure 2. Histograms of the quasar magnitudes (in Vega sys- 
tem) in all SDSS optical bands (ugriz) and UKIDSS near-IR 
bands(Y J H K) . Dashed lines indicate the median magnitudes. 



dian values of magnitudes and magnitude uncertainties are 
also indicated in Fig. 2 and Fig. 3. 

The parameters of this SDSS-UKIDSS sample of 8498 
quasars, including the coordinate, position offset, redshift 
and all magnitudes in SDSS and UKIDSS bands, are sum- 
marized in Table 1 (only a portion is shown here for dis- 
playing the content, and the whole table is available in the 
electronic version). 



3 THE COLOUR-Z RELATIONS 

With the SDSS-UKIDSS sample of quasars, we can calcu- 
late the median colour-z relations from u ~ g to H — K 
colours, which have important applications in determining 
the quasar selection criteria and photometric redshifts. In 
order to obtain reliable colours in both optical and near-IR 
bands, we only use the quasars with magnitude uncertainties 
in all 9 bands (ugrizY J H K) less than O.lmag. This reduces 
the quasar number from 8498 to 7400. Because the redshift 
range for our quasar sample is from 0.1 to 5.2, we are only be 
able to derive the colour-z relations for z < 5.2 quasars. We 
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Figure 3. Histograms of the quasar magnitude uncertainties 
in all SDSS optical bands (ugriz) and UKIDSS near-IR bands 
(YJHK). Dashed hnes indicate the median values of the magni- 
tude uncertainties. 



divide the 7400 quasars into different redshift bins with a bin 
size of 0.05, 0.1 and 0.2 for 2 < 3, 3 < z < 4.5 and z > 4.5 
respectively, and derive the median value of each colour at 
each redshift bin. We notice that the median u — g and g~r 
colour values derived from our sample are quite uncertain 
when the redshift is larger than 3 and 4.3 respectively, due 
to the small number (< 5) of quasars in the redshift bin. 
Therefore, we only give the median colour values for u — g 
at 2 < 3 and g ~ r at z < 4.3. 

In a previous work, Hewett et al. (2006) presented three 
tables of optical/near-IR colours derived for redshifted blue, 
average and red model quasar spectra, which are based on 
the quasar model SED given in Maddox & Hewett (2006). 
Chiu et al. (2007) plotted the colour-z relations generated 
from a modified version of the above quasar SED and 
compare them with the observational data of 2837 SDSS- 
UKIDSS quasars. In Fig. 4 we show the median colour- 
z relations derived from 7400 SDSS-UKIDSS quasars and 
the comparisons with the colours of all 8498 quasars in our 
present SDSS-UKIDSS sample. The colour-z relations given 
by Hewett et al. (2006) using the mean quasar spectrum 
are also shown for comparison. The agreement between our 
median colour-z relations and those in Hewett et al. (2006) 
is generally very well, especially for all colours in optical 
bands and Y — J colour in near-IR band. However, the dif- 
ferences of these two colour-z relations can be clearly seen at 
z < 0.5. Comparing with our median colour-z relations and 
the data, the colour-z relations of Hewett et al. (2006) seem 
to be redder in near-IR colours (except Y — J) and bluer in 
optical colours for z < 0.5 quasars. This is probably due to 
the continuum shape assumption in Maddox & Hewett et 
al. (2006) because the host galaxy contribution is significant 
for z < 0.5 quasars. Chiu et al. (2007) modified this con- 



tinuum shape and obtained the colour-z relation consistent 
with observations. 

In Table 2 we present our median colour-z relations. For 
u — g colour at 3 < z < 3.5, g — r colour at 4.3 < z < 5.0 and 
all other colours at 5 < 2 < 6.4, the values from the Table 
26 in Hewett et al. (2006) are also added for the purpose of 
determining the quasar selection criteria and photometric 
redshifts in the following sections. 



4 SEPARATING QUASARS AND STARS IN 
THE COLOUR-COLOUR DIAGRAM 

Separating quasars from stars is an important but difficult 
task in any quasar survey. Usually this can be done by the 
optical colour selection methods for low redshift quasars 
with z <2.3, as did by both 2DF and SDSS. At higher red- 
shifts, quasars usually have similar optical colours as stars 
so the contamination from stars to the quasar candidate se- 
lection is serious. However, as suggested by Warren et al. 
(2000), Chiu et al. (2007) and Maddox et al. (2008), this 
can be improved by incorporating the near-IR colours of 
quasars, because quasars are relatively more luminous in K- 
band than stars even though their optical colours are almost 
the same in some cases (e.g. for quasars with Z'^2.7). 

In some previous works, the colour-colour diagrams for 
quasars and stars were analyzed by including the UKIDSS 
near-IR colours. Hewett et al. (2006) proposed to separate 
quasars with 6 < z < 7.2 from dwarf stars in the i — Y 
vs. Y ~ J diagram. Maddox et al. (2008) studied the selec- 
tion criterion of 2 < 4 quasar candidates in the g — J vs. 
J — K diagram. Chiu et al. (2007) investigated the different 
colour-colour diagrams in optical and near-IR bands with a 
sample of 2837 SDSS-UKIDSS quasars, and found g — r vs. 
u — g diagram and H ~ K vs. J ~ H diagram are more effec- 
tive in separating quasars and stars than other colour-colour 
diagrams. They also proposed to use the Y ~ K vs. u ~ z 
diagram to select low redshift (2 < 3) quasars. 

In section 2, we have obtained a sample of 8498 quasars, 
which forms the largest sample of SDSS-UKIDSS quasars 
currently available. In order to study the quasar-star sepa- 
ration in the colour-colour diagrams, we also cross-identified 
8996 stars with both SDSS DR7 and UKIDSS/LAS DR3 
data following the same procedure as for our quasar sam- 
ple. These stars are classified as 'STAR' or 'STARJLATE' 
in SDSS spectroscopy. All magnitudes of stars in SDSS ugriz 
bands are converted to Vega magnitudes and corrected for 
the Galactic extinction. 

We have produced many colour-colour diagrams with 
our SDSS-UKIDSS quasar and star samples, including five 
colour-colour diagrams mentioned in the previous works 
mentioned above. Obviously separating quasars from stars 
is generally difficult in most colour-colour diagrams, but 
the colours involving K-band magnitude, such as H ~ K, 
,J — K and Y — K, can certainly provide some helps. The 
best diagram we found to separate quasars with 2 < 4 
and stars is the Y — K vs. g — z diagram, which is plot- 
ted in Fig. 5. From it, we see that a criterion defined by 
Y — K > 0.46 * [g — z) + 0.53, can well separate the 2 < 4 
quasars and stars. This criterion is obtained by searching 
for the highest efficiency in separating known SDSS-UKIDSS 
quasars and stars in our samples. With this criterion, we can 
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recover 8301 quasars in our sample of 8498 quasars (with a 
percentage of 97.7%), including 8277 of 8397 z < A quasars 
(with a completeness of 98.6%) and 24 of 101 z > 4 quasars 
(with a percentage of 23.8%). For 8996 SDSS-UKIDSS stars, 
only 204 of them (with a percentage of 2.3%) satisfy this cri- 
terion. Although the percentage of mis-identifying stars as 
quasars is very low, we must bear in mind that the actual 
contamination rate of stars could be higher than our esti- 
mation because SDSS only observed a tiny fraction of stars 
with particular colours spectroscopically. However, we be- 
lieve that our selection criterion is still fine as most SDSS 
untargetted stars (both early type and late type) should be 



within the stellar locus of our Y — K vs. g — z diagram and 
would not be selected as 2: < 4 quasar candidate by SDSS. 
In Fig. 5 we also show the predicted colours at different 
redshifts from the median colour-z relations of quasars (see 
section 3), which also confirm that our criterion is applicable 
for identifying quasars with z < 4. We notice that Chiu et 
al. (2006) has proposed a selection criterion in the Y — K 
vs. u — z diagram to select quasars, but their criterion can 
only be applied to 2: < 3 quasars because the Lya line moves 
out of the u-band for z > 3 quasars. In comparison, our new 
criterion has the advantage to select quasars with redshift 
up to 4. In addition, when we plot our SDSS-UKIDSS data 
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Figure 5. The Y — K vs. g — z colour-colour diagram of the SDSS- 
UKIDSS quasar and star samples. Green crosses, blue crosses and 
red circles denote quasars with z ^ 2.2, 2.2 < 2: ^ 4 and 2 > 4, 
respectively. Black dots and cyan dots denote stars classified as 
'STAR' and 'STARXATE' by SDSS spectroscopy Dashed line 
indicates the separation criterion Y — k > 0.46*(g— z)-|-0.53. Solid 
curve is derived from the median colour-z relation (see section 3) 
of quasars. The solid triangles, from left to right, mark the colours 
of quasars with z=0.02, 4, and 5 respectively, derived from the 
colour-z relation. Typical error bars are shown in the lower-right 
corner. 



Figure 6. The J — K vs. i — Y colour-colour diagram of the 
SDSS-UKIDSS quasar and star samples. Open triangles, open 
squares and solid squares denote M, L and T dwarfs given in 
Tables 9 and 10 of Hewett et al. (2006). Other symbols have the 
same meaning as in Fig. 5. Dashed line indicates the separation 
criterion J — K = 0.45 * [i — Y) + 0.64. Solid curve is derived 
from the colour-z relation (see section 3) of quasars. The solid 
triangles, from left to right, mark the colours of quasars with 
z=0.02, 5.3, 5.7 and 6.4 respectively, derived from the colour-z 
relation. Typical error bars are shown in the lower-right corner. 



of quasars and stars, we find that they are more separated 
in the Y — K vs. g — z plot than in the Y — K ys. u — z plot. 

Because the SDSS 2-band photometry is not as deeper 
as in ugri bands, our selection criterion in the Y — K ys. 
g — z diagram may be problematic for objects with z mag- 
nitude fainter than 20.5 (in AB system). Therefore, we also 
propose another criterion in the Y — K ys. g — Y diagram, 
namely, Y — K > 0.32 * — y) + 0.67, to separate quasars 
from stars. With this criterion, we can recover 8315 quasars 
in our sample of 8498 quasars, with a percentage of 97.8%, 
including 8276 of 8397 z < A quasars (with a completeness 
of 98.6%) and 39 of 101 2 > 4 quasars (with a percentage 
of 38.6%), but we also select 304 of 8996 stars (with a per- 
centage of 3.4%). The star contamination is higher than the 
case of using the criterion in the Y — K ys. g — z diagram. 
In the following study, we will mainly use the criterion in 
the Y — K YS. g — z diagram for selecting z < A quasars in 
order to keep a lower contamination rate of stars. 

For quasars with z > 4, obviously we need other selec- 
tion criteria. In most colour-colour diagrams, these high-z 
quasars have similar colours as late-type stars because the 
strong Lya line of quasars moves to the red part of optical 
spectra. By looking carefully at the various colour-colour 
diagrams, we realize that the high-z quasars and late-type 
stars often have different J — K and i — Y colours. We 
make sl J — K vs. i — Y colour-colour diagram (shown in 
Fig. 6) with the data of SDSS-UKIDSS quasars and stars. 
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Figure 7. Fraction of the 704 FIRST radio-selected 
SDSS/UKIDSS quasars recovered by applying two quasar 
selection criteria. The short and long dashed lines denote the 
cases of using the criteria in the Y — K/g — z and J — K/i — Y 
plots, respectively. The redshift bin size is 0.3. 
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Figure 8. The Y — J vs. i — z colour-colour diagram of the 
SDSS-UKIDSS quasar and star samples. Symbols have the same 
meaning as in Fig. 7. Solid curves are derived from the colour- 
z relation (see section 3) of quasars. The solid triangles, from 
left to right, mark the colours of quasars with z=0.02, 5.6, 6.0 
and 6.4 respectively, derived from the colour-z relation. Straight 
dashed lines indicate the selection criterion i — z > lA Y — J < 
0.1 * {i — z) + 0.25 for z > 5.6 quasars. Typical error bars are 
shown in the lower-right corner. 

Although some z < 4 quasars have the same colours as 
stars, the z > 4 quasars can be well separated from stars 
with a criterion J - K > 0.45 * (i -Y) + 0.64. In Fig. 6 
we also show the predicted colours from the median colour- 
z relations of quasars (see section 3), which also confirm 
that our criterion is applicable for identifying quasars with 
z < 5.3, although some quasars with 2.5 < z < 3, whose 
colours {i — Y ~0.5 and J — K ~0.8) are similar as stars in 
this colour-colour diagram, could be missed by applying this 
criterion. We check this criterion with our SDSS-UKIDSS 
sample of quasars and stars, and find that using the crite- 
rion J - K > 0.45 * (i - y) -I- 0.64 we can recover 8289 of 
8498 quasars in our sample (with a percentage of 97.5%), 
including 100 of 101 2 > 4 quasars (with a percentage of 
99.0%) and 8189 of 8397 z < 4 quasars (with a percentage 
of 97.5%). Only one z > 4 quasars and 208 z < 4 quasars 
are missed. For 8996 stars, 474 of them (with a percent- 
age of 5.3%) are mis-identified as quasars by our criterion. 
Therefore, the criterion in the J — K vs. i ~ Y diagram is 
better than the criterion in the Y — K vs. g — z diagram for 
identifying high-z quasars (with 4 < z < 5.3), although the 
latter is better for identifying z < 4 quasars. 

Perhaps a better way is to combine these two criteria 
for identifying quasars with redshift up to 5.3, namely, 

Y -K> 0.46 *{g~z) + 0.53, or (1) 

J - A' > 0.45 * (i - y) + 0.64. (2) 

With the above criteria, we can recover 8447 of 8498 quasars 
in our sample (with a percentage of 99.4%), including 8347 



of 8397 z < 4 quasars (with a percentage of 99.4%) and 100 
of 101 4 < z < 5.3 quasars (with a percentage of 99.0%). 
However, we would also mis-identify 566 stars as quasars, 
with a percentage of 6.3% of the total number of stars in our 
current SDSS-UKIDSS sample. In this case, we can achieve 
the highest efficiency in recovering both z < 4 and 4 < z < 
5.3 quasars, but the contamination rate of stars is also higher 
than using the individual criterion we proposed above. 

Because the radio quasars identified in SDSS were not 
selected by the colour-based methods (Richards et al. 2002), 
we can make a check with the FIRST radio-detected SDSS 
quasars to estimate the completeness of applying our selec- 
tion criteria, similar as what has been done in Richards et 
al. (2006). In our sample of 8498 SDSS/UKIDSS quasars, 
704 of them were detected by FIRST. With their photomet- 
ric data in SDSS and UKIDSS, we can check how many of 
them can be recovered by our quasar selection criteria in 
the Y — K/g — z and J — K/i ~ Y colour-colour diagrams. 
In Fig. 7 we show the fraction of recovered quasars at each 
redshift bin for two criteria we proposed above for z < 4 
and 4 < z < 5.3 quasars respectively. It is clear that with 
our Y — K/g — z criterion we can achieve the fraction higher 
than 95% in each redshift bin with z < 3.5, while with our 
J — K/i — Y criterion we can recover all FIRST quasars with 
3.5 < z < 5. Especially with our Y — K/g — 2 criterion, we 
can recover the quasars with 2 ~ 2.7 very well, which seems 
to be difficult in the case of using the SDSS quasar selection 
criteria alone (Richards et al. 2006). Such a check supports 
that with our two criteria we can identify quasars at a com- 
pleteness higher than 95% if their redshifts are smaller than 
5. 

Although our current SDSS-UKIDSS quasar sample 
does not include any quasar with 2 > 5.3, we can explore 
how to improve the high-z quasar selection by using the 
colour-z relations. For quasars with 2 > 5.3, especially those 
with 2 > 5.7, from Fig. 6 we can see clearly that the con- 
taminations of M, L and T dwarfs are much serious in the 
J — .K" vs. i — Y diagram. This is a problem already known 
for identifying high-z quasars (Fan et al. 2001b; Hewett et 
al. 2006), and the methods to solve it all involved near-IR 
colours. Fan et al. (2001b) proposed a selection criterion for 
2 ~ 6 quasars in the 2 — J vs. i — z diagram by requesting 
i — 2 > 2.2 and z ~ J < 1.5. Hewett et al. (2006) mentioned 
that 2 > 6 quasars can be identified in the i — Y vs. V — J di- 
agram as they have bluer Y — J colours than L and T dwarfs. 
Here we present another colour-colour diagram (Y — J vs. 
i — 2) in Fig. 8 and propose another selection criterion for 
high-z quasars: z — 2 > 1.4 and Y — J < 0.1 * (i ~ z) + 0.25. 
It is clear that using both i drop-out and near-IR colour can 
help us to separate 2 > 5.6 quasars from stars and low-z 
quasars, as well as T and L dwarfs effectively, though the 
contaminations from some stars are still unavoidable. Fur- 
ther spectroscopic observations on some high-z quasar can- 
didates selected with our criterion are expected. 



5 PHOTOMETRIC REDSHIFT ESTIMATION 

With the colour-z relations derived in section 3, we can es- 
timate the photometric shifts of quasars by comparing the 
observed colours with the predicted colours at certain red- 
shift by the colour-z relations. Combining SDSS ugriz and 
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Figure 9. Comparison of piiotometric redsiiift estimation results 
of 8498 SDSS-UKIDSS quasars for the cases of using SDSS colours 
only (upper panels) and using SDSS-UKIDSS colours (lower pan- 
els), and Zsp denote photometric and spectroscopic redshifts 
respectively. 



UKIDSS YJHK magnitudes, we can obtain 8 colours from 
optical u — g colour to near-IR H — K colour for the quasar 
candidates, then using the minimization method we can 
obtain the most probable photometric redshifts for them. 
Here the is defined as (see Wu et al. 2004): 



[(mi,, 



z ) {rrii, observed l^j, observed^ 



.(3) 



where the sum is obtained for all 8 SDSS-UKIDSS colours, 
mi,cz — rrij^cz is the colour in the colour-z relations, 
rrii^observed — Trij ^observed IS the observcd colour of a quasar, 

(Jmi absarvad ^'1'^ """"j obaerued ^''^ UUCertaiutieS Of Ob- 

served magnitudes in two SDSS-UKIDSS bands. 

We notice that the Lya emission line moves out of the u 
band and g band when the redshift is larger than 3.5 and 5 
respectively, so in our derived colour-z relation in section 3, 
we don't have the u~ g value for z < 3.5 and both the u~ g 
and g — r values for 3.5 < z < 5. This means that we can only 
use 7 and 6 of 8 SDSS-UKIDSS colours if the input redshift 
is larger than 3.5 and 5 respectively. In order to compare the 
values for the cases of using different number of colours at 
different redshifts, we actually use the 'x^ jN (where N is the 
number of colours and is 8,7,6 respectively for input redshift 
z < 3.5,3.5 < z < 5 and 2; > 5) instead of to determine 
the photometric redshift by obtaining the minimum olx^ /N 
at certain redshift. An IDL program is made for searching 
the photo-z of quasars by taking the above procedures into 
account. 

We estimate the photometric redshifts of 8498 SDSS- 
UKIDSS quasars with our program. In order to see more 
clearly whether using the UKIDSS colours can improve the 



photo-z estimation, we also estimate the photometric red- 
shifts of these quasars using the 4 SDSS colours only. Sim- 
ilarly as we did above, actually only 3 and 2 of 4 SDSS 
colours can be used for input redshifts larger than 3.5 and 
5, and photometric redshift is determined by searching the 
minimum of x? 1^- In Fig-9 we show the comparison of 
the photo-z results obtained with SDSS and SDSS-UKIDSS 
colours respectively. It is obvious that using both SDSS and 
UKIDSS colours can improve the photo-z accuracy substan- 
tially. For 8498 quasars, 7223 of them (with a percentage 
of 85.0%) have photometric redshifts given by the SDSS- 
UKIDSS colours consistent with the spectroscopic redshifts 
within lA^I ^ 0.2. 6102 of them (with a percentage of 71.8%) 
have |Az;| ^ 0.1. In comparison, if using the SDSS colours 
only, 6179 of quasars (with a percentage of 72.7%) have 
photometric redshifts consistent with the spectroscopic red- 
shifts within \^z\ ^ 0.2, and 4898 of them (with a per- 
centage of 57.6%) |Az;| 0.1. Therefore, by adding the 
UKIDSS YJHK data, the accuracy of photo-z can be im- 
proved with a percentage of 12.3% for \/S.z\ ^ 0.2 and 14.2% 
for |A2:| ^ 0.1, respectively. This substantial improvement 
is really what we expected. From Fig. 9 we notice that there 
are still large failures in the photo-z estimations for some 
quasars, even if the UKIDSS data are added. This is mainly 
due to the photo-z degeneracy (the X^/N minimum values 
are comparable at two or more different input redshifts) (Wu 
et al. 2004), and different colours of some special quasars, 
such as the reddened quasars (Weinstein et al. 2004). More 
dedicated approaches, such as using the additional morphol- 
ogy and colour information, could help to further improve 
the photo-z accuracy. 



6 TESTS WITH A SMALL UKIDSS/EDR 
SAMPLE AND VARIABLE SOURCES IN 
SDSS STRIPE 82 REGION 

In the two sections above, we presented the selection crite- 
ria of quasar candidates at different redshift and the pho- 
tometric redshift estimation method. Obviously more tests 
on these techniques need be done with some existing photo- 
metric and spectroscopic samples before their applications 
to a large quasar survey. 

First we test our selection criteria with a sample of 209 
objects classified as QSOs and 252 objects classified as A, 
K, and M stars in Maddox et al. (2008). These objects (with 
K4: 17.0) were selected from UKIDSS/LAS Early Data Re- 
lease and observed with the AAOmega multi-object spectro- 
graph. In Fig. 10 we plot these objects in the Y — K ys. g — z 
diagram and compare them with our quasar selection crite- 
rion, Y-K > 0A6*{g-z)+0.53. Among 209 quasars, 186 of 
them have z ^ 2.2 and 182 can be found with our criterion, 
23 of them have 2.2 < 2 < 4 and 19 of them meet our cri- 
terion. Therefore, using our proposed selection criterion we 
can recover 201 of 209 quasars (with a percentage of 96.2%) 
in Maddox et al. (2008). However, 13 of 252 stars, including 
one A stars, nine K stars and three M stars, also satisfy our 
criterion. The fraction of these mis-identified stars is 5.2% 
of the total number of stars and 6.1% of the total number 
of objects satisfying our selection criterion. Because all 209 
quasars in this sample have z < 4, this test confirms that our 
criterion is indeed applicable for identifying z < 4 quasars. 
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Figure 10. Testing the quasar selection criterion in Y — K vs. g—z 
colour-colour diagram with the observational data in Maddox et 
al. (2008). Red, black, and cyan triangles denote the A, K, and 
M stars respectively. Green and blue crosses denote quasars with 
2 ^ 2.2 and 2.2 < 2 ^ 4. Dashed line represents our quasar 
selection criterion shown in Fig. 6. Typical error bars are shown 
in the lower-right corner. 
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Figure 11. Comparison of the photometric redshifts with the 
spectroscopic redshifts for 209 quasars (left panels) in Maddox 
et al. (2008) and 819 SDSS identified quasars (right panels) in a 
sample of variable sources in Stripe 82. 



Figure 12. Upper panel: Y — K vs. g — z colour-colour diagram 
of all 4671 variable sources in SDSS Stripe 82 (Sesar et al. 2007) 
with SDSS-UKIDSS data. Green solid circles denote the 169 RR 
Lyrae star candidates. Lower panel: Only the variable sources 
with SDSS spectroscopy are shown. Green and blue crosses denote 
the quasars with z ^ 2.2 and 2.2 < 2 < 4 respectively. Black 
and cyan triangles denote the sources classified as STAR and 
STARXATE by SDSS respectively. Dashed lines in both panels 
represent the quasar selection criterion shown in Fig. 6. 



Maddox et al. (2008) indicated that 17 quasars in their sam- 
ple do not have SDSS spectra, therefore, this test also implies 
that with our criterion we are able to recover these SDSS 
missing quasars. In the left panels of Fig. 11 we show the 
comparison of the photometric redshifts estimated by our 
program (see section 5) with the spectroscopic redshifts for 
209 quasars in this sample. 158 quasars have |Az| ^ 0.2 and 
131 quasars have |A2;| ^ 0.1, with a percentage of 75.6% 
and 63.2% of the total 209 quasars, respectively. This also 
confirms the robustness of our photometric redshift estima- 
tions. 

The second test we made is for a sample of variable 
sources in SDSS Stripe 82, which is a region of 290deg^ lo- 
cated along the celestial equator in the southern Galactic cap 
and was repeatedly observed by SDSS photometry (Ivezic et 
al. 2007). Sesar et al. (2007) presented 13051 candidate vari- 
able sources with criteria g < 20.5, a{g, r) ^ 0.05mag and 
X^{g,r) ^ 3 (here a{g,r) and X'^(Si^) the variability 
amplitude and per degree of freedom in g and r bands 
respectively) from 1.4 million unresolved sources in Stripe 
82 observed at least 4 times in each of SDSS ugri bands. Fol- 
lowing the same procedure as in section 2, we cross-identify 
their near-IR counterparts with the UKIDSS/LAS DR3 and 
find 4671 of 13051 sources with detections in both Y and K 
bands. In the upper panel of Fig. 12 we plot the colours of 
these 4671 variable sources in the Y — K vs. g — z diagram 
and compare them with our quasar selection criterion. The 
SDSS magnitudes of these 4671 variable sources are taken to 
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Figure 13. Histogram of photometric redshifts of 3539 quasar 
candidates selected with our criterion from variable sources with 
g < 20.5 in Stripe 82. Red dashed histogram shows the spectro- 
scopic redshift distribution of 819 SDSS identified quasars. Green 
dotted histogram shows the redshift distribution of quasar can- 
didates with ivega < 19.1, which is the SDSS spectroscopic limit 
for z < 2.3 quasars. 



be the mean magnitudes given in Sesar et al. (2007). It can 
be clearly seen that these variable sources are divided into 
two categories approximately by our criterion. 3834 of these 
4671 sources (with a percentage of 82.1%) satisfy our crite- 
rion, while 837 sources, which do not satisfy our criterion, 
obviously locate in the stellar locus as we see in Fig. 5 and 
Fig. 10. In Fig. 12 we also plot 169 sources in this sample, 
which meet the selection criteria for RR Lyrae stars (Ivezic 
et al. 2005; Sesar et al. 2007). 145 sources of these 169 RR 
Lyrae candidates (with a percentage of 85.8%) locate in the 
stellar locus, and we believe that they are most likely real 
RR Lyrae stars. However, 24 of these RR Lyrae candidates 
meet our quasar selection criterion. Future spectroscopic ob- 
servations on these 24 RR Lyrae candidates will help us to 
clarify their nature. 

Among 4671 variable sources with both SDSS and 
UKIDSS Y and K band data, 834 of them have been 
identified by SDSS spectroscopy as quasars with redshifts 
0.1 < 2 < 3.3. In the lower panel of Fig. 12 we also compare 
their colours with our selection criterion. We can clearly see 
that 826 of them (with a percentage of 99.0%) satisfy our cri- 
terion, and only 8 quasars are exceptions (5 of them actually 
locate very close to our criterion line in Fig. 12). 72 of 4671 
sources have been identified as 'STAR' and 'STAR_LATE' 
by SDSS spectroscopy respectively. Two of them (with a 
percentage of 2.8%) satisfy our quasar selection criterion, 
and the rest 70 stars clearly locate in the stellar locus in 
Fig. 12. This again proves the robustness of our quasar se- 
lection criterion. Among 834 SDSS identified quasars, 819 
of them have UKIDSS data in all YJHK bands. We es- 



timate their photometric redshifts with our photo-z pro- 
gram and compare their photometric redshifts with spec- 
troscopic ones. The results are shown in the right panels of 
Fig. 11. 717 quasars have \Az\ ^ 0.2 and 627 quasars have 
lA^I ^ 0.1, with a percentage of 87.5% and 76.6% of the 
total 819 quasars, respectively. This tight agreement demon- 
strates again that using both SDSS and UKIDSS colours we 
can obtain more accurate photometric redshifts for quasars. 

Among 3834 quasar candidates selected with our selec- 
tion criterion in the Y — K vs. g — z diagram from 4671 
variable sources in Stripe 82, 3539 sources have all UKIDSS 
YJHK detections. We estimate their photometric redshifts 
with our photo-z program and the results are shown in Fig. 
13. 3519 quasar candidates have Zph ^ 3.5, 20 candidates 
have 3.5 < Zph < 5, and none candidate has Zph > 5. This is 
understandable because our selection criterion in the Y — K 
vs. g ~ z diagram is mainly for selecting z < 4 quasars (we 
also used our second criterion in the J — K vs. i — Y diagram 
to select possible z > A quasar candidates, but this only con- 
tributed 61 additional candidates). In Fig. 13 we also show 
the spectroscopic redshift distribution of 819 SDSS identi- 
fied quasars with all UKIDSS YJHK detections, and the 
photometric redshift distribution of our quasar candidates 
with i magnitude in Vega system iab < 19.1, which is the 
SDSS spectroscopic limit for z < 2.3 quasars. Comparing 
these two redshift distributions, we realize that even at the 
same magnitude limit, with our criterion we can probably 
select more quasars with z ^ 2 than SDSS. Especially for 
z ~ 2.7 quasars, which are largely missing in SDSS because 
of their similar colours as stars, can be at least partially 
recovered with our quasar selection criterion. For our 3834 
quasar candidates (with Y and K detections), only 826 of 
them were identified as quasars by SDSS spectroscopy, which 
is only a small fraction (21.5%) of our quasar candidates 
selected from variable sources in Stripe 82. This is mainly 
because the SDSS spectroscopic limit is iab < 19.1 for low- 
z quasars, which is shallow than gAB < 20.5 for selecting 
variable sources in Stripe 82 (Ivezic et al. 2007; Sesar et al. 
2007). We also notice that 207 sources of our quasar candi- 
dates in Stripe 82 have been recently identified as 199 new 
z < 3 quasars, 2 white dwarfs and 6 unclassified objects by 
the 2dF-SDSS LRG and QSO survey (2SLAQ) (Groom et al. 
2009). Therefore, deeper spectroscopy than SDSS, such as 
LAMOST quasar survey which is expected to reach a limit 
of lAB = 20.5, will hopefully produce a much larger sample 
of quasars in the future. 

In addition, we found that 837 of 4671 variable sources 
with both SDSS and UKIDSS (V and K) data do not satisfy 
our quasar selection criterion in the Y ~ K vs. g ~ z dia- 
gram. They clearly locate in the stellar locus in this diagram 
and most of them are probably stars. These 837 star-like 
sources only take a small fraction (17.9%) of our 4671 SDSS- 
UKIDSS variable sources in Stripe 82, which confirms that 
the variability study is very helpful for discovering quasars 
(Sesar et al. 2007). Our results suggest that with the near- 
IR photometric data we can further improve the accuracy of 
selecting quasars from optical variability study. This is prob- 
ably important for the variability study in the future large 
photometric surveys, such as that with the Pan-STARRS 
(Kaiser et al. 2002) and Large Synoptic Survey Telescope 
(Tyson 2002). 
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7 DISCUSSIONS 

Quasar candidate selections largely rely on the multi-band 
photometric colours. We have demonstrated that combin- 
ing the UKIDSS YJHK data with SDSS ugriz data wiU 
lead to significant improvements in both quasar selection ef- 
ficiency and photometric redshift accuracy. With our quasar 
selection criteria involving both SDSS and UKIDSS colours, 
we can probably select more z > 2 quasar candidates than 
SDSS even at the same spectroscopic magnitude limit. We 
believe that at least some SDSS missing quasars with z ~ 2.7 
can be recovered with our selection criteria. Deeper spectro- 
scopic observations on the quasar candidates selected with 
our criteria are expected to confirm our results and produce 
a more complete and much larger quasar sample in the fu- 
ture. 

The quasar selection criteria we proposed in this pa- 
per are fairly simple, involving only cuts in the two-colour 
diagrams for optical point sources. More complicated ap- 
proaches, such as the multi-dimensional search, would be 
more optimal. To check this, we have tried several differ- 
ent ways to search the criteria for separating quasars and 
stars, including 3-dimensional colour-colour diagrams, sup- 
port vector machine (SVM) and Bayesian classification al- 
gorithm, but we did not find significant improvements. The 
detailed comparisons of the results by adopting these differ- 
ent approaches with our current work will be presented in a 
separated paper. 

The study presented in this paper is mainly based on the 
SDSS-UKIDSS sample of 8498 quasars. Although the cross- 
identifications were done by finding the closest UKIDSS 
counterparts of the SDSS quasars and the offsets between 
the SDSS and UKIDSS positions are mostly within 0.5", the 
mis-identifications of some sources are unavoidable. How- 
ever, by checking the colour-z relations and colour-colour 
diagrams of these SDSS-UKIDSS quasars, we believe that 
the fraction of such mis-identifications should be very low. 
Therefore, the main results obtained in this paper will not 
be affected by the mis-identifications of some sources. 

Because there is no quasar with z > 5.3 in our SDSS- 
UKIDSS quasar sample, our quasar selection criteria can 
only be applied to identify z < 5.3 quasars. Although we 
present a criterion in the Y — J vs i ~ z diagram for finding 
z > 5.6 quasars based on the colour-z relation given by 
Hewett et al. (2006), this criterion obviously needs to be 
tested with a large sample of z > 5.6 quasars with both 
SDSS and UKIDSS data. Moreover, for the high-z quasars 
with 2 > 5, the photometric redshift estimation will largely 
rely on the UKIDSS colours because both u — g and g ~ r 
colours are not available. We also need a sample of jz > 5 
quasars with both SDSS and UKIDSS photometric data to 
check the accuracy of our photometric redshift estimations. 
This is expected to be done in the future when a larger 
high-z quasar sample is available. 

Although we made two tests using a UKIDSS/EDR 
sample and a sample of variable sources in SDSS Stripe 
82 to check the robustness of our quasar selection criteria 
and the photometric redshift method, obviously we still need 
more tests. We have selected some quasar candidates with 
Jab < 19.1 using our criteria from the variable sources in 
Stripe 82, which do not have SDSS spectroscopy, and will 
identify them with optical spectroscopy. This will provide 



direct evidence for whether using our SDSS/UKIDSS selec- 
tion criteria can help to find missing quasars in SDSS. On 
the other hand, finding quasars with certain criteria and 
comparing them with those in the existing spectroscopic 
surveys may not be enough for testing the completeness, as 
most known quasars were identified with the colour selection 
methods. A better way to do it may be to run simulated 
colours of quasars, with assumptions on their continuum, 
emission line and reddening, for different redshift and lumi- 
nosity (Fan 1999), and check how many simulated quasars 
can be found with the selection criteria. This will be studied 
in our future work. 

The Chinese LAMOST quasar survey will be done in 
the next a few years with the most efficient spectroscopic 
telescope in the world, and is expected to obtain a much 
larger quasar sample (with the magnitude limit of i=20.5 
or 21) than the previous surveys. The results of our current 
study, if confirmed with further spectroscopic observations, 
will be helpful for preparing the input catalog of quasar can- 
didates and estimating their photometric redshifts for the 
LAMOST quasar survey, especially in the overlapped sky 
area between SDSS and UKIDSS. By doing necessary tests 
of various selection criteria with the spectroscopic observa- 
tions in the pilot survey, we will determine the best crite- 
ria for quasar candidates selections and apply them to the 
LAMOST quasar survey. Although the UKIDSS/LAS will 
cover only a sky area of 4000deg^ , with the near-IR colours 
we will be able to select more complete quasar samples with 
redshift up to 5 than previous surveys. This large and com- 
plete quasar sample will be very important to many fur- 
ther studies such as those on the quasar luminosity function, 
quasar clustering, and large scale structure in the universe. 
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Table 2. Median colour-z relations of SDSS-UKIDSS quasars 
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Quasar candidate selection and photometric redshift estimation 



Table 3. Table 2 continued 
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Note: All colours arc derived from magnitudes in Vega system. The u — g 



colour at 3 < 2 < 3.5, g — r colour at 4.3 < z < 5.0 and other colours at 
5 < 2 < 6.4 are adopted from Tabic 26 of Hewett et al. (2006). 
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